Optimal local expansion of W states using linear optics and Fock states 
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We derive the maximum success probability of circuits with passive hnear optics for post- 
selectively expanding an A''-photon W state to an (A'' + n)-photon W state, by accessing only one 
photon of the initial W state and adding n photons in a Fock state. We show that the maximum 
success probability is achieved by a polarization-dependent beamsplitter and n — 1 polarization- 
independent beamsplitters. 
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I. INTRODUCTION 

Multipartite entanglement has rich structure and di- 
verse properties, which can be exploited for many kinds 
of applications. With an increasing number of parties 
sharing entanglement, the fact that the structure be- 
comes more complex also means that preparation of such 
an entanglement may become more difficult, especially if 
one tries to generate the whole state in a single step. 
Instead, one may start with an initial entangled state 
among a small number of subsystems, and then expand 
it by adding ancillary subsystems. In fact, many in- 
teresting classes of multipartite entangled states, such 
as Greenberger-Horne-Zcilingcr (GHZ) states [1], clus- 
ter states [2] and W states [3], are defined for an ar- 
bitrary number of qubits, which makes the expansion 
strategy look attractive. There are several proposals of 
state expansion for multipartite entangled states by lo- 
cal manipulation on a single site without accessing other 
qubits [4-6]. For example, in the case of GHZ states, 
a deterministic local expansion of an A^-qubit state to 
an [N -\- n)-qubit state is possible in principle. In op- 
tical systems, several experimental demonstrations have 
been performed [7-10] by using quantum parity checking 
gates [11]. 

W states are an interesting class of multipartite en- 
tangled states in that they have a web-like entangle- 
ment structure. An A'^-qubit W state is represented by 
IWat) = (|10---0)-H|010---0) + --- + |0---01))/V77. Ev- 
ery qubit in [Wat) has bonds with every other qubit, and 
the pairwise entanglement survives even if all the other 
{N — 2)-qubits arc discarded [3, 12, 13]. In recent years, 
there have been a number of theoretical proposals and ex- 
perimental demonstrations using W states in multiparty 
protocols such as quantum key distribution [14], leader 
election [15] as well as preparation of W states in optical 
systems [16-26]. In the case of W states, a deterministic 
local expansion is impossible even in principle because 
the marginal state of the remaining untouched N — 1 
qubits is different for [Wat) and |Wiv+n), so it is worth 
seeking nontrivial ways of expanding W states probabilis- 
tically. Since the qualitative difference from GHZ states 
in the expandability may arise from the difference in the 
nature of multipartite entanglement, study of efficient lo- 
cal expansions of W states is interesting theoretically as 



well as practically. 

For photonic polarization-based qubits, recent propos- 
als and a demonstration include the expansion of [Wat) 
to |Wjv+i) [4, 5] and to |Wjv+2) [6, 27]. These expansion 
schemes are composed of passive linear optics and one 
or two photons in a Fock state. In this paper, we ad- 
dress the question of what is the best way of expanding 
[Wat) to \W N+n) for photonic polarization-based qubits. 
We discuss the maximum success probability of local ex- 
pansion methods composed of passive linear optics and 
an ancilla mode in an n-photon Fock state. We derive 
the maximum success probability, and also show that 
it is achieved by a polarization-dependent beamsplitter 
(PDBS) and n — \ polarization-independent beamsplit- 
ters (BSs). In the case of n = 2, the optimal success 
probability is higher than that of the expanding gate pro- 
posed so far [6]. 

This paper is organized as follows: In Sec. II, we de- 
scribe the expansion schemes considered in this paper. In 
Sec. HI, we derive an upper bound on the success prob- 
ability for expanding W states. In Sec. IV, we explicitly 
construct a linear optical circuit that achieves the de- 
rived upper bound on the success probability. Finally, in 
Sec. V, we give a brief summary and conclusions. 



II. EXPANSION METHODS OF W STATES 

The optical circuit considered in this paper for expand- 
ing W states is composed of only passive linear optics, 
that is, of PDBSs, BSs, phase shifters, wave plates, and 
arbitrary linear losses. As shown in Fig. 1 (a), the circuit 
has two input spacial modes 1 and 2. One photon from 
an A'^-photon polarization-entangled W state, which is 
represented by \Wn) = (|VH • • • H) + |HVH • • • H) + • • • + 
|H • • • HY))/y/N where |H(V)) is the state of an horizon- 
tally (H-) and vertically (V-) polarized photon, is fed to 
mode 1, and an H-polarized n-photon Fock state (we de- 
note it by [tiH)) is fed to mode 2. We define the successful 
operation of the circuit to be the events where the pho- 
tons come out from the output spacial modes 1, . . . , n+1, 
one by one, namely, exactly one photon from every out- 
put mode. We require that the output state is exactly 
the {N + n)-photon W state for the success events. 

As is well known, linear optical losses can be equiva- 
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FIG. 1: (Color online) Local expansion circuits of jWjv) to 
|Wjv+n)- (a) The expansion circuit considered in this paper. 
This circuit is composed of only passive linear optics, (b) The 
expansion circuit in Fig. la is equivalent to a lossless linear 
optical circuit like this figure. Solid(dashed) arrows mean 
H(V)-polarized modes. 



lently described by lossless optical circuits vifith auxiliary 
spacial modes. This is because a linear loss with trans- 
mission T is equivalently realized by using a BS with 
transmission T. Since a BS is a lossless component with 
two input modes and two output modes, we can simulate 
losses with lossless components by introducing auxiliary 
input modes initially in the vacuum and the same num- 
ber of auxiliary output modes from which the 'lost' pho- 
tons escape. The circuit in Fig. 1 (a) is thus described 
equivalently by a lossless circuit with L{> n + 1) input 
and output spacial modes. The family of schemes con- 
sidered in this paper is defined as follows. The expan- 
sion circuit has 2L input and output modes composed 
of L spatial modes with each having two polarizing (H- 
and V-polarized) modes. One photon from an iV-photon 
polarization-entangled W state enters into input spacial 
mode 1 and an H-polarized n-photon Fock state enters 
into input spacial mode 2. All of the other input spacial 
modes receive vacuum states. Arbitrary lossless linear 
optical operations are applied to the 2L input modes. 
We show this expansion circuit in Fig. 1 (b). 



Let Ivac) 



i(out) 



be the state of all the 2L input(output) 



modes in the vacuum. The unitary operation U for 
any passive linear optical operations satisfies U\vac) in = 



I uac) out, and relevant actions of U are represented by 

L 



Ua\^U^ - E(^^h61h + /5jv6]v) , (1) 

L 

Ua\yU^ = ^(7,h6]h+7,v6Jv) , (2) 



and 



(3) 



where ajjjj.y^ and (i!2Yi{y) H(V)-polarized pho- 

ton creation operators for input spacial mode 1 and 
2, ^jH(v) ^^'^ H(V)-polarizcd photon creation opera- 
tor for output spacial mode j, and ajH(v)j Pjii(v) ^^nd 
7jH(v) E^re complex numbers. Using the commutation re- 
lation [ajs,a],^,] = [bjs,bys'] = Sjj>5ss'{s,s' = H,V) and 
Eqs. (1) - (3), we obtain 



(4) 



for Q, n' = a, P, 7. By denoting [Hat) = |H • • • H), the in- 
put state to the expansion circuit in Fig. 1 can be written 
as 



(5) 



+ |HAr_i) 

while state N+n) which we desire as an output is 

|WAr+„ 



I ( "+1 

= Vn~1\Wn-i) ® n 



n+1 n+1 

b]^\vac)out 

1=1 j^i 



vac)ont 



(6) 



The post-selected events where the photons come out 
from output spacial mode 1, . . . ,n + 1 one by one are 
described by a projector written as 



n 



post — X! ' ■ ■ X! 

si=H,V s„+i=H,V 

/n+1 \ /n+1 \ 

I ) K'«c>(«ac|out f n^''vj -(7) 



^/=l 



By applying U and Hpost to Eq. (5), wc should obtain 
Eq. (6). Hence 



npostC/|WA.)|nH) /f^|WAr+„), 



(8) 
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where Pguc is the probabihty of success. By operating 
(Wjv-il ® (i;ac|out(nr=i'^/H) on Eq. (8) from the left, 
we have 



, nUN + ii) 

VPsuc = \l Vo , 



(9) 



where 



VO ^ (vacUl^Ubmy^^^^^^lvac),^ (10) 



n+l ri+1 

/3jH Yl "iH 
1=1 j^i 



Here we have used Eqs. (l)-(3), (5) and (6). Similarly, 
by operating (Hat^iI ® {vac\out{hv ^(h) on Eq. (8) 



from the left, we have 



n\{N + n) 



N 



where 



n+l 



Tji = {vac\ont biv Y\_bi}i\u 



'^1v('^2h)^ 

n! 



(12) 



|wac)i„(13) 



n+l 



n+l n+l 



= 7»v Yi ^jti + '^»v X! ^i'^ n "'^H ' 

for 1 < I < n + 1. Here we have used Eqs. (l)-(3), (5) 
and (6). From Eqs. (9) and (12), we obtain 



P., 



n\{N + n) I 



\m\ 



and 



Vo 



Vn+l 



(15) 
(16) 



III. BOUND ON THE SUCCESS PROBABILITY 

In this section, we derive an upper bound on success 
probability Pguc- Since the cases with Pguc = are irrel- 
evant for the upper bound, we focus on the cases with 
^suc > here. The requirement of producing the exact 
W state, Eq. (8), narrows down the choices of the unitary 
operator U significantly. As shown in Appendix A, in or- 
der to satisfy Pguc > 0, it is necessary that a^v = Piv = 
for 1 < i < 71 + 1 and Jlj^^i '^jH 7^ 0. Thus, from 
Eqs. (11) and (14), we are allowed to express r/o and 
Vi, for later use, as 



^n+l 



n+l 



.0 = (^n«.HjE- 

n+l 

Vi = 7jV Y[ "jH 



(17) 
(18) 



for 1 < i < n + 1. Define Pi = |q;,hP (1 < « < n + 1), 

P = {Pi: ■ ■ ■ , Pn+l), 



n+l 



and 



5(P)^5]P, , 



n+l 

H(P) EE Yl P, 

J = l 



(19) 



(20) 



(11) Because n"=i "jH and S{P) < 1 from Eq. (4), we 
obtain 

PeTZ = {P\P, > 0(1 < i < n + 1), < S{P) < 1}. (21) 



Combining a^v = ftv ~ for I < i < n + 1 
with EjLila^n/^jH + a*vl3jv) = and Ej=i(l"iHp + 

l«.vP) - Et=i(lftHp + = 1 from Eq. (4), we 

have 



n+l r, 
3 = 1 •' 



= E E (o.-a*.("+m-. 

j = l s=H,V 

< E E \Qs^a;,{n + l)\' 



3 = 1 s=H.V 

n+l 



E^-(" + i)^ 



(22) 



where Cm = ^^^^ CjV = for 1 < i < n + 1, and 
Oh = CiV = for n + 2 < i < L. Here we have used the 
Cauchy-Schwarz inequality. From Eqs. (17), (20) and 
(22), we obtain a bound on |?/oP written as 



<P(P)^n(P) (^E]^ + 



(23) 



Note that F{P) > since we are focusing on 
the cases with Pguc > 0. Using Eqs. (16), (18), 



(19) and E?=il7jvp < 1 from Eq. (4), we have 



\vo\'s{p) = E]ZiP3\v3\' = E;:iP3\i3y\'m;-Pk = 

E"=i l7jVpn(P) < n(P). Then we obtain another 
bound on bnP written as 



\vo\'<G{P)^ 



_H(P) 



S{P) 



From Eqs. (15), (21), (23) and (24), we obtain 

Psuc < — — max H(P) , 
N Pen ^ ' ' 



(24) 



(25) 



where H{P) = min{P(P), G(P)}. 

Before conducting the optimization over P in Eq. (25), 
let us discuss physical intuition behind the bound H{P). 
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Recall that Pi = laiHp is the probability of a photon in 
the input mode 2 in Fig. 1 (a) to appear at the output 
mode i. One of the bounds on |?7oPi C!{P), was derived 
through constraints on , which is proportional to the 
probability of having a V-polarized photon at the output 
mode i and n H-polarized photons in the other n output 
modes, one in each [See Eq. (13)]. From the definition of 
G{P) in Eq. (24), we see that 



and 



SiP) = Pi 



Pi 



Pn 



Pn 



+ 1 



(26) 



(27) 



give the maximum of G{P). This means lossless and 
equal distribution of the n H-polarized photons incident 
on the input mode 2 is the best for maximizing the am- 
plitude of the terms including one V-polarized photon. 
This result does not change even if we regard photons as 
classical distinguishable particles, since the origin of the 
V photon (input mode 1) and that of an H photon (in- 
put mode 2) are uniquely determined and no interference 
occurs. 

On the other hand, the bound F{P) stems directly 
from a constraint on |?7oPi which is proportional to the 
probability of (n -I- 1) H-polarized photons to appear at 
the {n + 1) output modes, one by one. In this case, there 
are (n + 1) indistinguishable paths, depending on which 
of the {n + l) output photons is traced back to the input 
photon in mode 1. As a result, the total amplitude 770 is 
given by the sum over {n + 1) terms as in Eq. (11), which 
can be rewritten as 



^70 



(28) 



Let us see the interference among these terms at the 
choice of P satisfying Eqs. (26) and (27). The param- 
eters f3iii{i = 1,...,L) describe how the input photon 
in mode 1 is distributed. They must satisfy the unitar- 
ity condition of Eq. (4), which gives, under the lossless 
condition of Eq. (26), a constraint 



n+l 



(29) 



Together with Eq. (27), we see that the total amplitude 
770 always vanishes regardless of the choice of the param- 
eters /3iii{i = l,...,n-t-l). This is the reason why we 
have F{P) ^ for P ^ {{n + 1)-^, ...,{n + 1)-^). In- 
cidentally, the case with n = 1 is equivalent to the well- 
known two-photon interference effect at the symmetric 
beamsplitter [28], in which the two photons never leave 
separated. The above result with general n can thus be 
regarded as an extension of the (1, l)-photon case to the 
(l,n)-photon case. What is interesting here is that the 
symmetry is required only for the n-photon input, and 



not for the one-photon input, to achieve the complete 
destructive interference. 

In order to obtain a nonzero value H{P) > 0, one must 
go away from the point Pi = • • • = P^+i = {n + 1)~^ by 
dropping either the lossless condition Eq. (26) and/or 
the symmetry condition Eq. (27). If one breaks the sym- 
metry, the orthogonality condition Eq. (29) no longer im- 
plies 770 = in Eq. (28). On the other hand, introduction 
of loss relaxes the orthogonality condition Eq. (29) itself, 
since it is equivalent to introducing auxiliary modes in 
Fig. 1 (b), namely, L > n + l. The condition Eq. (29) 
then changes to 



71+1 



L 



L 

E 

i=n+2 



iYPiV 7 



(30) 



which allows more freedom in the choice of parameters 
/3iH(i = l,...,n + l). 

In the following, we optimize over P in Eq. (25) by 
deriving necessary conditions for P to achieve the maxi- 
mum of H{P). We assume Pi < • • • < Pn+i without loss 
of generality, and we consider the two cases, S{P) < 1 
and S{P) ~ 1, separately. 

In the case of S{P) < 1, for P to be a local maxi- 
mum, there exists e > such that H{P + APo) < H{P) 
for any APq = auQ + l3vo with + (3'^ < e, where 
Uq = (1,0,..., 0) and Vq = (0, . . . , 0, 1). As shown in 
Appendix B, this leads to 



{uo,VF){vo,VG) = {vo,VF){uo,VG) 



(31) 



where {X,Y) means the inner product between X and 
Y. From V,P(P) = dF{P)/dP^ = p-\F{P) - 
U{P)Pr^),dG{P)/dPi = G{P){P-^ - S{P)-^) and 
Eq. (31), we obtain 



1 



P 



n+l 




J2pA = 0- (32) 



J=2 



Since Pi < ••• < P„+i, Eq. (32) means Pi = ••• = 
P„+i S{P){n + 1)^1. At this point, P(P) and G(P) 
are regarded as functions of single parameter S{P). We 
show in Appendix C that S{P) = 1 — {71 + gives the 
local maximum of H{P), whose value is 



lossy 



_ n"(?7 + 2)" 
{n + l)3"+i 



(33) 



In the case of S{P) = 1, since Pi > for all i, for P 
to be a local maximum under the constraint S{P) = 1, 
there exists e > such that H(P + APi) < H{P) 
for any APi = aui + /3vi with + < e, where 
Ml = Mo - (Mo,V5(P))(V5(P),V5(P))-iV5(P) = 
{n{n + l)^^,-{n + l)"\...,-(n + 1)"^) and Vi = 
vo - {vo, VSiP)){VSiP), V5(P))-iV5(P) = (-(?! + 



1)-^ , . . . , -{n + l)-\n{n + 1)-^). This leads to 
{ui,VF){vi,VG) = {vi,VF){ui,VG) 



(34) 
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from Appendix B. From Eq. (34), wc obtain 



n+l 



n+l 



Here the first factor is nonzero since Pi < Pn+i from 
F{P) > 0. In Eq. (35), since {P^^ - Pf^) > and 
(P,7_^_\ - P~^) < for all j, we obtain either P^ = Pi or 
Pj = Pn+i for every j, which implies Pi — ■ ■ ■ = Pm < 
Pm+i = ■ ■ ■ = P„+i(l < m < n). We thus find that, at 
the local maximum, 

F{P) = pr~'p:^r[^Pn+i 

+ (n + 1 - m)Pi - [n + ifPiP^+i] , (36) 



and 



with 



G{P) = P^Pl'.+l-"' 



mPi + (n + 1 — m)Pn+i = 1 



(37) 



(38) 



In Appendix D, we derive the local maximum of H{P) 
under the constraint S{P) ~ 1, whose value is i?m in 
Eq. (D8). Furthermore, in Appendix E, we show the 
maximum of Hm is given by m = 1, and its value is 



where 



PI 



,opt _ 2n + 3 - V4n + 1 



2(n2 + 2n + 2) 



(39) 



(40) 



From Eqs. (33) and (39), maxp^-jz H{P) is equal to 
max{7Jiossy, -ffi}- In Appendix E, we show Hi > i?iossy, 
and thus we obtain 

Ps..<^^^('-Yprii~prr. m 
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FIG. 2: (Color online) A linear optical circuit which achieves 
the upper bound in Eq. (41). 




FIG. 3: The dependence of the maximum success probability 
fVnax in Eq. (42) on n and A^. 



one by one with equal probability, that is, they are set as 
Tfc = (?i-fc)(n+l-fc)-i and Rk = (n+l-fc)~^ [29]. In 
this case, am = \/7h, a2H = • • • = a(,i+i)H = y/n-^Rn, 



ajuin + 1) and 7jv = ctjn for all j. All 



of the other variables of q;jH(V)5 /^jH(V) a-nd 7jH(V) a-re 
equal to zero. We have rjQ — ■ ■ ■ — in this case, 

and the success probability of this circuit is equal to the 
righthand side of Eq. (41). Hence we conclude that the 
maximum of Pour is 



Prr 



n\{N + n) (1 
N 



P°'"^(l — p°P*)" 



(42) 



IV. EXPLICIT CONSTRUCTION OF AN 
OPTIMAL CIRCUIT 

Here we construct an optical circuit which achieves the 
righthand side of Eq. (41). The circuit in Fig. 2 is com- 
posed of a PDBS and n~l BSs, and this circuit has n+l 
input and output spacial modes with no auxiliary spacial 
modes. We post-select the events where exactly one pho- 
ton comes out from every output mode. We denote the 
transmittance and reflectance for H(V)-polarized pho- 
tons of the PDBS by Pkcv) and Ph(V)5 and the trans- 
mittance and reflectance of the fc-th BS (BSfc) by Tk and 
i?fc(l < fc < n — 1). The parameters of the PDBS are set 

as Th = Pv = P7\ Ph = Tv = 1 - ^T'- The param- 
eters of the BSs are set to output n photons from n BSs 



The dependence of Pmax on n and N is shown in Fig. 3. 

In the case of n = 1, we see Pmax = [N + l)/(5iV) for 
Th(v) = (5 ± a/5)/10 and Tv(h) = (5 T a/5)/10. This 
value is the same as that in Ref. [4, 5]. In the case of 
n 2, we have P^ax = + 2)/(125A^) for Th = 4/5, 
Tv = 1/5, Ti = 1/2 and Pi = 1/2. This value is higher 
than [N + 2)/(16A^) given by the circuit composed of 
only half BSs in Ref. [6]. 

We also construct optical circuits which achieve the 
success probability of the other local maximums, Pm in 
Eq. (D8) of Appendix D and Piossy in Eq- (33). A circuit 
which achieves Psuc = n!(A^ + n)N~^H„i is constructed 
with a PDBS and n - 1 BSs. {BSi, BS„i_i} and 
{BS,„, . . ., BS„+i} are placed at each output of the PDBS 
in series, as shown in Fig. 4 (a). This circuit has n+l 
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FIG. 4: (Color online) Linear optical circuits which achieve 
(a) P,uc = n\{N + n)N-^H„„ and (b) P^uc = n!(iV + 



input and output spacial modes without auxiliary spacial 
modes. Parameters of the PDBS and the BSs are set as 

Th = Rv = fn^m, Rh = Hy ^ 1~ m^,„ where .Jm is in 
Eq. (D4), Tk = {m-k){m + l-k)-^ and Rk = {m + 1- 
k)-^ for 1 < A: < m - 1, and Tk = {n - k){n + 1 - k)-^ 
and Rk — (n + 1 — fc)^^ for m < /c < n — 1. In this 

case, am = • • • = q;,„h = V'n^^Tn, a(m+i)H = • • • = 

-1 



and 



ajH(ri + 1) 



7jV 



a^H for all All of the other variables of 



Q^jH(v)j /3jH(v) and 7jH(v) are equal to zero. 

A circuit which achieves Psuc = flossy = nl{N + 
n)N~^H\ossy is constructed with a PDBS and n BSs. The 
PDBS and BSi, . . ., BS„ arc placed in scries, as shown in 
Fig. 4 (b). We post-select the events where the photons 
come out from output spacial modes of n BSs, one by one. 
This circuit can be regarded as a lossless circuit with n+2 
input and output modes, but if we regard it as a circuit 
with 2 inputs and n+l output, as in Fig. 1 (a), it is a lossy 
circuit. Parameters of the PDBS and the BSs are set as 
i?,H = 1 - (n + Th = (n + i?v = 0,Tv = 1, 

Tfc = {n + l- k){n + 2- k)-^ and Rk ^ {n + 2 - k)~^. 
In this case, am = ■ ■ ■ = Q(„^i-)h 



a(n+2)H 

I < j < n 
7jv = V(n 



= y'RH{n+l)-\ 

R}i{aJ^ - aj}i{n + I)) for 
-V^"(«+2)h("- + 1), and 
for 1 < j < n + 1. All of the other 



1 



■(n+2)H 



variables of ajH(v) > /^jH(V) and 7jH(V) are equal to zero. 

Let us compare the success probability flossy for the 
lossy circuit with the global maximum flnax for a few 
examples. In the case of n = 1, we see fliossy = 3(iV + 
l)/(16iV) for Th = 1/4 and fl = 1/2, whereas fl^^ax = 
{N + 1)/{5N). In the case of n = 2, we have Piossy = 
128(iV + 2)/(2187Af) - 0.059(^ + 2)/N for Th = 1/9, 
Ti = 2/3, T2 = 1/2, while P^ax = 0.064(iV + 2)/N. 



V. CONCLUSION 

We have derived the maximum success probability of 
the circuits composed of passive linear optics and an an- 
cilla mode in an n-photon Fock state for post-selectively 
expanding an iV-photon polarization-entangled W state 
to an ( -f ri)-photon polarization-entangled W state, by 
accessing only one photon of the initial W state. Whereas 
the symmetry in W states suggests that photons from 
beamsplitters should be equally distributed among n + l 
output modes, bosonic nature of photons requires us to 
introduce either optical losses or to break symmetry in 
order to reduce interference effects between one photon 
from the A^-photon W state and the photons from the 
Fock state. In fact, both cases possess local maximums 
at which the success probability does not increase by 
infinitesimal changes in variables. We showed that the 
overall maximum success probability is achieved by a 
PDBS and n—1 BSs. In the case of n = 2, the maximum 
success probability is higher than that of the expanding 
gate proposed in Ref . [6] . 
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Appendix A 



In the following, we prove that when Pguc > is satis- 
fied, wc have a^v = and /3iv = for 1 < i < ri -|~ 1, and 
ttjH 7^ 0. For nonzero success probability Psuc > 0, 
it is necessary to have 770, f?i 7^ 0. From Eqs. (5) and (8), 
the term in IWat) including must not be transformed 
into the terms including bly by the unitary operator U. 
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Then we obtain 



T.., ^ {vacU I n C/-"ihI«2h) 



n+1 



n 

n+1 



ajH + a<v X! ]J a/cH = , (Al) 



for 1 < i < n + 1 and 

'n+1 



(Wac out 



n u 




^ /3,v n "fcv = . (A2) 

3 = 1 k^j 



IS 



Because J2i=i ^i^^m Hi^i "(v = ™]o n"=i ajv = 
obtained from Eqs. (Al) and (A2), we have 0^=1 '^jv — 
0. Assuming that aiy = 0, we find three facts: (i) 
n"^/"jH 7^ because 77/ = livlY-^l ^ 0, (ii) 
= for 1 < j < n+1 because rji = jiy YVj^l aja 7^ 

and 



n+1 



n+1 



vac) in 



for 1 < I) < n+1 from Eq. (8), and (iii) /3(v = from 
r; = Av n",!/ "iH = 0. Therefore, from iTj^i "jv = 
and recursive use of (i)-(iii), we have aiy = and Piy ~ 
for 1 < i < n + 1, and n?=i 0!jH 7^ 0. 



Appendix B 



and (a, VG) = (6, VF) > 0. Hence both (a + b, VF) > 
and (a + b, VG) > 0, which imphes the direction of 
a + b increases F and G at the same time. This fact 
contradicts with H{P + AP) < H{P) for any AP. In 
the case of {u,VF){v,VG) < {v,VF){u,VG), both 
(-a - b, VF) > and (-a - b, VG) > are satisfied, 
which also contradicts with H{P + AP) < H{P) for any 
AP. Therefore we obtain Eq. (Bl). 



Appendix C 

We derive the maximum of H{P) = min{F(P), G(P)} 
in the case of Pi = • • • = P„+i = S'(P)(n + 1)~^, namely, 



F{P) = , \! (I - S{P)) 



and 



G(P) = 



S{P) 
S{P) \n + l 



n+1 



(CI) 



(C2) 



F{P) and G{P) are regarded as functions of single pa- 
rameter S{P). From Eq. (CI), we have 



dF(P) (S{P) 



dS{P) 



n+1 



(n-5(P)(n+l)). (C3) 



We find that F{P) decreases monotonously for n(n + 
l)-i < 5(P) < 1 from Eq. (C3), and G(P) increases 
monotonously for < S{P) < 1 from (C2). From 
Eqs. (CI) and (C2), the solution of P(P) = G(P) is 
given by S{P) = 1 - (n + 1)"^ Since n{n + 1)"^ < 
l-{n + l)-2 is satisfied, S{P) = 1 - (n + 1) ^ gives the 
maximum of H(P), whose value is 



lossy 



n"(n + 2)" 
(n + l)3"+i 



(C4) 



Here we prove the following statement. Suppose 
H{P) = min{P(P), G(P)}. Let e > be a constant and 
u and V be arbitrary vectors. If Fl(P + AP) < H{P) 
for any AP = au + /3v with + < e, then 

(tt, VF) {v, VG) = {v, VF) {u, VG) , (Bl) 

is satisfied where {X, Y) means the inner product be- 
tween X and Y . 

We define two vectors a and b as 



and 



a = {u, VF)v - {v, VF)u , 



b= {v,VG)u- {u,VG)v . 



(B2) 



(B3) 



Suppose (m, VF)(t>, VG) > {v,VF){u,VG). From 
Eqs. (B2) and (B3), we obtain (a, VF) = (b, VG) = 



Appendix D 

We derive the maximum of H{P) = min{F(P), G(P)} 
under Eqs. (36) - (38), namely, 

F(P) = F(e) ^ r-'C-"'{mC + {n + l-m)^ 

-{n+l)'^C), (Dl) 
G(P) = G(e) = ^'"(^'i+i-m ^ (J32) 



with 



c 



1 — TTl^ 
n + 1 — 771 



(D3) 



Because ^ < < ^ < (rt + 1) ^ is satisfied. From 
Eqs. (Dl) and (D2), F(^) = G(0 has four roots, which 
are given by ^ = 0, tti"^ and the two roots of !{£,) = 
mC+{n+l-m)^^ {n + 1)2^C - CC = 0. Since /(O) > 
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and /((n + 1) < 0, there is only one root satisfying 
F{^) = G(0 for < C < (71 + which is 



2m,{n + 1) + 1 - y/Am{n + 1 - m) + 1 



2TO((n+ 1)2 + 1) 
From Eqs. (Dl), (D2) and (D3), wc obtain 



(D4) 



rn — 2/-n—m — l 



(C-0/(0 



where 

fiO 

and 



n + I — m 



n + 1 + 2nm 
n + 1 — m 



1 — m 



(D5) 



(D6) 



(D7) 



From Eq. (D7), we find that G(^) increases monotonously 
for < ^ < (ri + 1)^^. Since /(^) is a convex function, 
/(^m) = (e-l)(C-l)-l<Oand/((n + l)-i) = -2(n + 
1)^^ < assure that F(^) decreases monotonously for 
£.m < S, < {n + l)^^ . Therefore ^ = gives the maxi- 
mum of H{P), and the value is 



1 - m^n 

n+1 -1 



n+1— m 



(D8) 



Appendix E 



We show that Hi > -ffm(2 < m < n) and Hi > -ff lossy, 
where i?™ is defined in Eq. (D8) and iJiossy is in Eq. (33). 



Since H^ > 0, let us analyze the property of logi?m as 
a continuous function of m for 1 < m < n. Let h(m) = 
Hj^^dH,n/dm be its derivative. Because we have 



dh{m) 
dm 



(n + 1) - 2m 



m{n + 1 - m){Am[n + 1 - m) + 1)^/^ 



.(El) 



we find that m = (n + l)/2 gives the maximum of h(m). 
From d/i((n+ l)/2)/dn = 2(n + l)-i((?i+ 1)^ + l)-3/2 > 
for all n and h{{n + l)/2) — > {n od), we obtain 
h{m) < for all m. Thus, because dH„i/dm < holds, 
Hm takes its maximum for m = 1, whose value is 



Hi = {-] a(i-ar, 



(E2) 



from Eq. (D8), where is in Eq. (D4). 

Next we show Hi > i?iossy Since Hi > H(„+i)/2 is 
satisfied, wc show H(„+i)/2 > ^^lossy From Eqs. (33) 
and (D8), 



H(n+l)/2 



_ (n + l)4 
(n + 1)4-1 



(n + 1)4-1 (71 + 1)2-1 



> 1,(E3) 



is satisfied. Hence we obtain Hi > H] 



lossy • 
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